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£^ ' Abstract. We present a detailed comparison of three different methods designed 

to tackle nonequilibrium quantum transport, namely the functional renormalization 
group (fRG), the time-dependent density matrix renormalization group (tDMRG), 
and the iterative summation of real-time path integrals (ISPI). For the nonequilibrium 
single-impurity Anderson model (including a Zeeman term at the impurity site), we 
£f) ' demonstrate that the three methods are in quantitative agreement over a wide range of 

parameters at the particle-hole symmetric point as well as in the mixed-valence regime. 
We further compare these techniques with two quantum Monte Carlo approaches and 
the time-dependent numerical renormalization group method. 
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1. Introduction 



Nowadays possibilities of miniaturization down to the nanometer scale allow for studying 
the single electron transport in ultra-small devices such as, e.g., artificially designed 
molecular quantum dots [H El IH HJ [5] or nanotubes, with nontrivial emergent physics 
such as Fermi-liquid behaviour in quantum dots [6|. Quantum many-body systems 
driven to nonequilibrium tend to approach a stationary state being distinct in character 
from their ground state properties [TJ IB]. The details of the stationary state may 
in principle depend on the way how the system is driven out of equilibrium and 
on the character of the correlations in the system. There is a broad variety of 
interesting effects when a large bias voltage beyond the regime of linear response is 
applied to the system due to interaction effects or due to the nonequilibrium conditions 
0, HO, HI]- These features range from coherent transport properties, such as, e.g., 
resonant tunneling, to the Kondo effect. A prominent nontrivial model to study the 
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interplay of correlation effects and quantum transport is the single impurity Anderson 
model [12]. This model is believed to capture essential features of experiments with 
quantum dots [131 E]- Applying a finite bias voltage across the impurity region allows 
one to study nonequilibrium effects. This regime is defined by the breakdown of the 
standard approach of linear response theory when the transport voltage becomes too 
large. 

There are different approximative approaches to deal with the nonequilibrium 
situation. For example, quantum transport through a quantum dot in the Kondo 
regime has been described theoretically by Fermi-liquid theory for the asymptotic low- 
energy regime [15], by means of interpolative perturbation schemes [IE], by exploiting 
the integrability of the model [T71 US) , and in terms of the non-crossing approximation 
[l9j|20]. Moreover, there exists a large class of renormalization group approaches which 
are based on an expansion in the (renormalized) system-reservoir coupling and therefore 
allow for a treatment of the weak-tunneling regime. For example, the Anderson dot 
in the charge fluctuation regime [21] and the time evolution of the spin-boson model 
[22] have been initially considered in an early application of the real-time RG (RTRG) 
method in its real-time formulation [23]. Nonequilibrium extensions of the perturbative 
renormalization group [2U [25j [26] and of the flow equation approach [27] as well as the 
field-theoretical RG approach focusing on the Callan-Symanzik equation [28] have been 
applied to the nonequilibrium Kondo effect. A more refined analytical perturbative 
scheme based on a reformulation of the RTRG in frequency space (RTRG-FS) [29] 
provides a systematic and self-consistent treatment of observables [30] and correlation 
functions [31] in the nonequilibrium stationary state. Using RTRG-FS, approximate 
analytic solutions have been proposed for the nonequilibrium Kondo model for weak 
spin fluctuations [30| |3T| [32] and for the interacting resonant level model for weak 
or strong charge fluctuations [33J. Additionally, the RTRG approach also allows for 
a description of transient dynamics towards this state [22j [32]. A complementary 
treatment of transport properties through the Anderson dot can be based upon an 
expansion in the Coulomb interaction strength [3^1 35J. However, bare perturbation 
theory appears to be insufficient for larger values of the interaction, and therefore, it 
should be combined with the concept of the renormalization group. The latter prescribes 
how to renormalize interaction vertices in the nonequilibrium case. This is achieved 
by combining the Keldysh formalism with the functional renormalization group (fRG) 
approach in the quantum field theoretical formulation [361 E3 ESI [39] . Very recently, 
this scheme has been improved by accounting for the frequency dependence of the 
two-particle vertex [38] in analogy to corresponding developments of the equilibrium 
Matsubara fRG [401 W\\ . As a result of this improvement, the exponentially small scale 
of the Kondo temperature has partly been observed, namely in the nonequilibrium 
Fermi-liquid coefficient, that is in the second derivative of the self-energy with respect 
to bias voltage. Moreover, fRG has been applied to the study of transport properties of 
a Kondo quantum dot as well [42J. 

As an alternative to these approaches, numerical methods have been advanced 
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very recently to properly describe nonequilibrium quantum transport. One important 
technique is the real-time quantum Monte Carlo method (QMC) [43, 44J, which has 
been extended to non-equilibrium |45j HH1 H7J HEJ EH]. Although these results are 
valuable since they are numerically exact, the related calculations are limited to short 
to intermediate simulation times and the dynamical sign problem often does not allow 
one to reach the regime of the stationary current for low temperatures. Recent real-time 
QMC simulations [49] seem to indicate that the splitting of the Kondo resonance which 
was theoretically predicted earlier in (201 El] is an artifact of the perturbative approach, 
consistent with observations from tDMRG [50]. A novel technique based on quantum 
Monte Carlo (QMC) simulations with complex chemical potentials has been introduced 
in [51]. This is achieved by an analytic continuation to complex times by employing 
Matsubara voltages. Then, standard equilibrium quantum Monte Carlo techniques are 
used in this QMC variant [HTJ E2] • Beyond raising conceptually fundamental questions, 
a nonmonotoneous behaviour of the conductance on the bias voltage has been reported 
[5Tj . whose origin remains unclear and might be attributed to the (arbitrary) choice of 
the form of the scattering states on which the formalism relies. 

Moreover, the multi-layer multi-configuration time- dependent Hartree formalism 
has recently been introduced [531 |54], which is based on a clever decomposition of 
the overall Fock space and which has been applied to vibrationally coupled electron 
transport. In addition, the numerical renormalization group (NRG) method has been 
very successful in solving quantum impurity problems in equilibrium for several years, 
see |55j for a topical review. This approach has been generalized to nonequilibrium 
by using a discrete single-particle scattering basis [561 E0 EH1 Ell [601 EE] , resulting in 
the development of the time-dependent numerical renormalization group (TD-NRG). 
In addition, in [62] and for certain limits of the single impurity Anderson model, 
a dynamical steady state characterized by correlation induced current oscillations is 
presented by means of the time dependent density functional theory (TD-DFT). 

Furthermore, the adaptive time-dependent density matrix renormalization group 
method (tDMRG) that generally allows for the simulation of the real-time evolution of 
pure states in strongly correlated one-dimensional systems [631 El] has been applied 
to transport in nanostructures as well [SH ESI ESJ EZ] and most notably in the 
present context, to the calculation of current- voltage characteristics of quantum dots 
[501 E3 EHl HQ] . The steady-state currents are computed by taking suitable time-averages 
over time windows, in which initial transient effects have disappeared already. 

Finally, yet another numerical approach has been proposed that is based on the 
deterministic iterative summation of real-time path integrals (ISPI) [71]. It determines 
a Keldysh generating function for the time-dependent nonequilibrium current and yields 
numerically exact results since the systematic errors within the scheme are eliminated 
exactly. In a first step [71] , ISPI was developed for the single- impurity Anderson model 
and its validity has been confirmed by a detailed comparison with various approximative 
approaches within their range of applicability. 

It is the aim of this work to compare three of these methods, namely, tDMRG, fRG 
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and ISPI in detail, focusing on steady-state currents as a function of gate potential, 
voltage bias, temperature, and magnetic field. The goal is to establish their reliability 
by revealing an excellent quantitative agreement between results obtained with these 
techniques. We shall further check our results against those of other methods, namely 
the TD-NRG [59] and two QMC schemes, first the one from [511 152] and second, the 
real-time QMC from [ITJHS]. In section |2j we introduce the single-impurity Anderson 
model, which is one the most fundamental models for the description of quantum dots 
and which we will thus use to benchmark our methods. In section [31 we briefly review 
the fRG, tDMRG and ISPI schemes. The comparison of the respective numerical results 
and a discussion are contained in section HI 



In what follows we consider the single-impurity Anderson model p2] given by the 
Hamiltonian (h = 1) 

H = Hdot + Hi ea ds + Ht 



Here, Eq„ = Eq + aB with a =f, 4-= ± is the energy of a single electron with spin a on 
the dot. It can be varied by tuning a back gate voltage or by means of a Zeeman term 
oc B, under the assumption that the electron dispersion in the leads is not influenced 
by the magnetic field. The creation/annihilation operator for the dot electron is d\jd a , 
n a = d\d a has eigenvalues n a — 0, 1, and the on-dot interaction is U . The energies of the 
noninteracting electrons in the lead p = L/R = ± with chemical potential fi p = peV/2 
are denoted by e^ p (with fermionic operators c^ pa ). The dot is connected to the leads 
via the tunnel coupling t p . 

Within the fRG and the ISPI method, we assume the leads to be in thermal 
equilibrium at temperature T (kg = 1) and moreover take the standard wide-band 
limit with a constant density of states, p{€f), around the Fermi energy, yielding the 
hybridization T p = 7rp(ei?)|i p | 2 . In this work, we assume symmetric contacts, i.e., 

= Tr = T/2 and we scale all quantities with respect to T. The observable of 
interest is the symmetrized tunneling current I = (II — Ir)/2 with 



2. Model 



^ Eoafia- + Un^hi + 2j( 
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kpcr 




t , 



(2) 



kpa 



where I p {t) = —eN p (t) with N p (t) = (J2ka c k P a c kpa)t- The asymptotic long-time limit 
gives the stationary steady-state dc current as / = lim t ^ 00 /(t). 
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3. Methods 

Next, we briefly review the three techniques that we want to compare. 
3.1. fRG 

A detailed explanation on how we apply the fRG to the single-impurity Anderson model 
is given in reference [37]. Here we only sketch the main ideas. 

Spectral as well as transport properties of the model can be comprehensively 
described if an interaction contribution to the single-particle self-energy £ is somehow 
established [19J. A basic idea of the fRG approach consists of a formulation of a flow 
equation for S accounting for its full frequency dependence (therefore this RG approach 
is called functional) . The starting point of the flow is the noninteracting limit where 
the solution is known exactly. During the fRG flow the form of £ is continuously 
transformed, and the solution of an interacting problem is achieved when the flow 
terminates. 

Technically, this procedure is set up by artificially making the bare propagator g 
depend on a flow parameter A. Being functionals of the bare propagator, the interacting 
Green's and vertex functions acquire a dependence on A as well, which is described by an 
infinite hierarchy of coupled flow equations [72]. In this paper we focus on the flow of the 
one-particle irreducible vertex functions [73] I7i] which has proven to provide a successful 
approach to the physics of several low-dimensional correlated electron problems [7_5l |7_6] . 
The one-particle irreducible n-particle vertex function can be defined diagrammatically 
as the sum of all one-particle irreducible diagrams with n amputated incoming lines 
and n amputated outgoing lines. Here we refer to diagrammatic perturbation theory 
on the Keldysh contour which results from an expansion in powers of the two-particle 
interaction U. 

We truncate the hierarchy of flow equations by neglecting the flow of the three- 
particle vertex function. The flow equations for the self-energy E A and the two-particle 
vertex function 7 A then read 

— rSl/M = —~ 7l'9'l19 Snlr,/ (3) 



where indices occurring twice in a product imply summation over state and Keldysh 
indices and integration over independent frequencies. Furthermore, G x is the full single- 




A AAA ^ A AA AAA 

1'2'|34 ^3\3' ^4]4' 73'4'|12 + ^~~7l'4'|32 [^3|3' ^4|4' + ^3|3' ^4|4'] 73'2'| 
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(4) 



particle propagator and S x = —G x \d(g x )~ l / dX^ G x denotes the so-called single scale 
propagator. 

As flow parameter we use the hybridization constants T p by setting 
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where A flows from infinity to zero. This choice for the flow parameter ensures the 
validity of causal properties and of the fluctuation dissipation theorem in the special 
case of thermal equilibrium [77]. The latter is important in order to reproduce the Fermi 
liquid behaviour of the model near the particle- hole symmetric point. We approximate 
the coupling of the three channels in the flow equation (jlj) of the two-particle vertex by 
reducing the influence of each channel onto the other two to a renormalization of the 
interaction strength U. This leads to a frequency dependence of the two-particle vertex 
of the form 

7 A (n, X, A) ~ v + <^(n) + <£{X) + ^ A (A), (6) 

where II = oj\ + u 2 = oj[ + uj' 2 , X = u' 2 — to\ = uj 2 — to[, A = u[ — oj\ = oj 2 — oj 2 
are bosonic frequencies corresponding to energy exchange in different channels. In 
particular, II refers to the particle-particle and X to the exchange particle-hole channel, 
while A indicates the direct particle-hole channel. A similar approximation has been 
investigated in studies of the single impurity Anderson model based on the equilibrium 
Matsubara fRG |41j . In ([6]), v is the bare interaction vertex which is the initial value 
at the beginning of the flow, and ^123 are approximations to the corresponding parts 
of 7 produced by the flow equation, see [37]. This approximation considerably reduces 
the number of sampling points in frequency space which are required for the numerical 
solution of the flow equation. 

3.2. tDMRG 

The adaptive time-dependent DMRG method [631 164"] allows for the simulation of the 
time-evolution of pure states \ip(t)) at zero temperature (note that generalizations to 
mixed states and finite temperatures are possible, for a review, see [T8j and references 
therein). The basics of the technique itself are described in the original publications 
[63~ll64"] and the review [79]. In a nutshell, the key ingredient of the method is to represent 
the wave-function in a truncated, but optimized basis instead of the exponentially large 
basis of the full Hilbert space (see [78] for details). The application of tDMRG to 
transport in nanostructures was pioneered in [5Ql EEJ ETJ EHJ EH EQ] • We use the set-up 
of [501 EZ]) which we shall briefly summarize. 

In contrast to ISPI and fRG, we here apply tDMRG to a real-space representation 
of the Hamiltonian, i.e., we replace the terms in (CQ) describing the leads and the 
hybridization by 

ifleads + H T = -t p 2j (dl C l,l,a + h -C.) 
l=R,L;cr 

w Nl,r 

--4- ( C lna C l,n+la + h.C.) . (7) 

n=l;l=R,L;a 

Moreover, we work on finite systems with open boundary conditions. Num is the 
number of sites in the left (right) lead, with iV = Nl + Nr + 1. The hopping matrix 
element in the leads is W/4, where W is the bandwidth. In the present case of a 
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semielliptical density of states, we obtain T = 2tp/(W/4) for the hybridization. Note 
that the leads can also be modeled with other density of states by a different choice of 
(typically site-dependent) hopping matrix elements in the leads. A prominent example 
is the logarithmic discretization which is the standard choice in NRG, recently used in 
some tDMRG simulations as well [69, 80J. 

A simulation then starts by computing the ground state of ([7]), defining the initial 
state \ipo)i an d then, at time t = 0, the system is driven out of equilibrium by applying 
on-site energies ±eV/2 in the leads, mimicking the bias 



For the time-evolution \ip(t)) = exp(—iHt / h)\^ ) , we use a Trotter- Suzuki breakup of 
the time-evolution operator with a time-step of 5t = 0.1 in units of 1/(W/A). During 
the time-evolution, we measure the symmetrized tunneling current /(£), ([2]). Some 
typical results for the real-time data are shown in figure [TJ As a function of time, I(t) 
goes through a transient regime before a steady-state is reached. The transient time, 
besides its dependence on T (and, in the Kondo regime on the Kondo temperature 
Tjc), is proportional to 1/V [50] • The steady-state regime can be accessed in tDMRG 
simulations if the typical transient times are smaller than N/vp, where Vp is the Fermi 
velocity in the leads. At present, time scales of about 10/T can easily be reached at 
intermediate bias voltages as figure [U illustrates (see also the discussion in [50]). The 
final result of a simulation is, for our purpose, the steady-state current that is computed 
by taking time-averages over the real-time data. Methods to identify the quasi-steady- 
state and for finding the correct time-windows for averaging are discussed in [50l 166] 
and we refer the reader to [50], 1661 [81] for a discussion of the time scales governing the 
transient regime. 

The numerical errors inherent to the technique depend on two control parameters, 
first, the time-step St and second, the discarded weight Sp [78]. The latter is a measure 
of the truncation that the wave-function has been subjected to. In the example studied 
here, the discarded weight is the relevant quantity [50J. Moreover, finite-size effects 
need to be taken into account and thus runs are carried out for N < 96 at small biases 
V/t < 0.4 and N < 64 at large biases. All results presented in this work have been 
extrapolated in 1/N. 

Unlike ground-state DMRG [78], there is no criterion that favourably limits the 
entanglement encoded in the time-evolved wave-function (see also the discussion in [50] 
and references therein), rendering the simulation increasingly costly at long times and 
in particular, for large biases V > W/2. Qualitatively, the larger the entanglement, the 
larger the dimension of the truncated basis set needs to be that is used to approximate 
\ijj(t)). This is at present the most crucial limitation of the method. Nevertheless, 
tDMRG has been very successful in yielding current-voltage characteristics for two of 
the most basic models used to describe quantum dots, the interacting resonant level 
model [68] and the single-impurity Anderson model [50J. 




(8) 



n=l 



n=l 
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3.3. ISPI 

The ISPI scheme [71] is an in principle numerically exact method to deal with non- 
equilibrium transport through correlated quantum systems. The method is based on 
the evaluation of a Keldysh generating function [82] , which includes appropriate source 
terms to generate the observables of interest. Here, this is the nonequilibrium transport 
current. A real-time path integral 



Z[v) = Jv 



(9) 



with Grassmann fields (d, d, c, c) (here we use the same symbols for the Grassmann 
fields and the fermionic operators in ([I])) is constructed for the generating function. 
The action S corresponds to the Hamiltonian (CQ) and an external source term is added 
to the action that allows for computing the dc current at measurement time t m via 
I{tm) = — i-jr In Z[v] ■ The full real-time evolution is computed by successively 

^ ri=0 

inserting t m at different Trotter slices, e.g. in [71], figure 5 therein. Next, the path 
integral is Trotter discretized along the Keldysh contour. The interaction (U-) term is 
decoupled via a discrete Hubbard-Stratonovich transformation on each single Trotter 
slice, by which an auxiliary "spin" -field is introduced. For the single- impurity Anderson 
model, it can assume the values ±1. Then, all fermionic fields can be integrated out 
analytically. The remaining task is to perform the summation over all configurations of 
the (bosonic) auxiliary fields, reminiscent of a discrete path summation. Hence, in 

^W = EIlHdetG; 1 [M^]), (10) 

{s} a 

(with the total effective Keldysh Green's function G^ 1 ), the summation over all 
configurations s is performed numerically in a deterministic way. The key ingredient for 
this is the fact that the fermionic environment (the leads) induces non-local correlations 
in time which decay exponentially in the long-time limit for any finite temperature. 
This implies the existence of a characteristic memory time up to which all correlations 
are fully taken into account and beyond which they can be safely neglected for larger 
times, similar to the case of the spin-boson model [83]. This allows us to disentangle 
the complete path sum and to construct an iterative scheme to evaluate the generating 
function. By means of a properly chosen extrapolation procedure, both systematic 
numerical errors of the scheme, namely the Trotter error and the error due to the finite 
memory-time, are eliminated. If the extrapolation scheme will be convergent, one is 
left with the desired numerically exact values for the observable at a fixed measurement 
time. Convergence is generally problematic when both T and V approach zero. For 
more details, we refer the reader to |71j . 
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4. Results and Discussion 

For the comparison of the three approaches, we compute the symmetrized steady-state 
current (j2J). A comparison of the fRG and the tDMRG results has already been shown in 
[50] for zero temperature and both at particle-hole symmetry and in the mixed valence 
regime. Consistency between these two methods could be demonstrated for U/Y < 8 and 
eV > T K , where T K is the Kondo temperature. We now set out to include ISPI results 
into this comparison, starting with the zero temperature transport at the symmetric 
point E = —U/2 at finite on-dot interaction U. We shall then extend our study to the 
mixed valence regime Eq ^ —U/2, and end up with including a finite magnetic field and 
finite temperatures, presenting new data from all three methods. 

In figure El we show the results for the particle-hole symmetric point, E = —U/2, 
of the single-impurity Anderson model. For the case of U/Y = 2 we find an excellent 
agreement between the three methods and with the results from the nonequilibrium 
real-time QMC approach of [17]. The fRG is expected to be very reliable for U/Y = 2 
and we can thus use fRG to benchmark the other methods in this limit. Generally 
(the same applies to the following figures as well), the deviations at large bias voltages 
between tDMRG on the one hand and fRG and ISPI on the other hand are due to 
the finite bandwidth used in the tDMRG calculations, as compared to the wide-band 
limit taken in the other two techniques (see [50] for details). For the large bias voltage 
(around eV/Y ~ 4) tDMRG and ISPI deviate from each other by up to 4%. The high 
accuracy of the ISPI data is illustrated by our findings shown in figure [2} where the 
method agrees very well with both the fRG and the tDMRG, to be specific deviations 
are below 1.5%. For U/Y — 4 we find good agreement of the fRG and the tDMRG with 
the real-time QMC approach from [49]. In the case of U/Y = 8, fRG and tDMRG yield 
very similar results whereas the real-time QMC gives a slightly larger current. Yet these 
deviations are within the numerical error inherent to the techniques, see e.g. [50] . 

Next, we compare the results of fRG, ISPI and tDMRG at intermediate interactions 
and set U/Y = 3 (ISPI and fRG are computed at T = T/25.6). The respective current- 
voltage characteristics /(V) is shown in figure [3] (top). We find a nice agreement over 
a large range of V, while again, small deviations between ISPI and fRG occur at very 
large voltages. In figure |3] (bottom), we compare the differential conductance obtained 
from our three methods with the result of Han and Heary [51] (see figure 2 in that 
work). We numerically approximate dl/dV ~ [I[V + AV)] — I(V)/AV by choosing a 
sufficiently small AV. 

Strong differences between the QMC result and our three techniques occur and 
the non-monotonous features of the QMC data are not reproduced in our simulations. 
Moreover, the results of fRG and ISPI coincide at small to intermediate voltages, while 
some deviations occur at larger voltage 0. Within the error bars of the numerical 
differentiation and also taking into account uncertainties in the bare I{V), tDMRG 

| More recent versions of the method of [5T| have given somewhat better agreement [84 . An extensive 
discussion of this method has recently been presented in [55] . 
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yields results for dl/dV that are consistent with those from fRG and ISPI. 

In a next step, we further increase the Coulomb interaction and choose U/T = 5, 
for which TD-NRG results are available [59]. Using ISPI in its present formulation, this 
choice for U/T does not yet allow us to obtain converged results. In figure IU we show 
the results for / from fRG and tDMRG, together with the TD-NRG data by Anders 
[59] . We find an excellent agreement between fRG and tDMRG for eV > 1.75I\ The 
TD-NRG data are very similar to fRG at small bias voltages, while TD-NRG gives larger 
values for the steady-state current at large values of V than the other two methods. The 
origin of this discrepancy is, at present, not understood, yet we have reason to speculate 
that discretization errors may play a role here. tDMRG simulations for an Anderson 
impurity coupled to Wilson leads, i.e., applied to the same Hamiltonian that NRG 
treats, also result in too large currents at eV > Tk [69]. Using improved discretization 
schemes might provide a remedy for this problem and should be tested in future work 
(see, e.g., [Ml EH EH]). 

Next, we address the mixed valence regime for the case of U/T = 2. In figure El we 
show the steady-state current I for different single electron energies E ^ —U/2. The 
results from fRG and ISPI match perfectly from small bias voltages eV/T ~ 0.2 up to 
the strong non-equilibrium regime. Note that by construction, it becomes increasingly 
cumbersome (and finally impossible) to obtain converged ISPI results in the limit of 
vanishing bias voltages and low temperatures [71] . since then the correlations do not 
decay sufficiently fast enough to be truncated. In the present set-up, tDMRG has a 
tendency to overestimate the currents in the mixed valence regime, for more details, see 
[50J. This, we believe is the reason for slight deviations from the other two methods 
(see eV ~ 1.5r, see figure [SI lower left panel). Overall, the agreement between the three 
methods is obviously still very good, even away from the symmetric point. 

Moreover, we consider the case of a finite magnetic field B at the symmetric point, 
see figure EJ We find an excellent agreement between fRG and ISPI over the full range 
of bias voltage (deviations are below 1%) and the tDMRG data agree well with the 
two other methods for eV < 1.5r, where we can mimic the wide band limit by keeping 
r, U, B, eV < W. 

Finally, in the case of finite temperatures T (see figure [7]), for which at present 
no tDMRG is available yet, we again find a remarkable agreement between fRG and 
ISPI from the equilibrium up to the strong non-equilibrium case, with deviations below 
1.5%. Only at large temperatures both methods tend to disagree slightly for larger bias 
voltages and deviate up to 10%. We note that for increasing bias voltages, the fRG 
approach becomes more and more reliable and even asymptotically exact as V — > oo. 
ISPI also profits from growing bias voltages as they cut off the time non-local correlations 
in the leads' Green's functions more and more efficiently 

In the remaining part, we wish to point out advantages and disadvantages of the 
three techniques this work was mainly concerned with, the fRG, tDMRG, and ISPI. 
fRG is a computationally rather cheap tool to compute current-voltage characteristics 
over the entire bias range and at practically all temperatures. The approximations 
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taken in this method have been estimated to be valid up to U/T ~ 6. In particular, 
in [38] it has been shown that they yield a quantitatively accurate description of 
the linear conductance as a function of gate voltage, temperature and magnetic field; 
good descriptions of equilibrium properties of the model have already been observed in 
Matsubara implementations of the fRG [JH EE] • Concerning nonequilibrium properties, 
we obtain very good agreement for the nonlinear current and differential conductance 
for U/T < 8 with tDMRG data close to the most critical point E - U/2 = T = B = 
[50] . For sufficiently low (but still finite) temperatures, ISPI yields results that also 
compare well with fRG. Within the Keldysh fRG framework, we are able to recover 
|38] the nonequilibrium Fermi-liquid relation [15] and identify an exponentially small 
energy scale of the Kondo temperature in the corresponding Fermi-liquid coefficient 
which is the second derivative of the self-energy with respect to bias voltage (although, 
true Kondo physics for large U/T cannot be described in general [3E])- However, the 
reliability of fRG results at larger values of U/T > 6 obtained in the scheme truncated 
at the level of the three-particle vertex is unknown at present, since the justification to 
truncate the fRG hierarchy of equations can only be given by perturbative arguments. 
In the single-impurity Anderson model, it is not a priori clear that the contribution of 
a three-particle vertex is negligible at large U. This constitutes an interesting subject 
for future research. 

Using tDMRG, one can access the finite voltage range with eV > Tk, as has 
been demonstrated for U/T < 8 in [50] and in this work, while the full current-voltage 
characteristics can be obtained at small U/T ~ 4. The main limitation of the method is 
the entanglement growth, rendering long time-scales difficult to access at large voltages, 
which at present excludes an analysis of the deep Kondo regime. A partial solution 
to this problem is the use of leads with a logarithmic discretization, and then, the 
correct result is obtained for the linear conductance and U/T < 7 [69]. tDMRG results 
can be rendered quite accurate as well, yet since the method requires an exponentially 
increasing computational time both for reaching long time scales and to obtain very 
accurate data, in practice, one has to accept a finite numerical error, which can be 
estimated by tuning control parameters. The greatest advantage of tDMRG is that 
existing codes can be directly applied to both complex interacting structures as well as 
interacting leads [90] . 

Regarding the ISPI scheme, an important advantage is that whenever it converges 
with respect to its internal control parameters, numerically exact results are obtained, 
where no sign problem restricts the accessible simulation times and the steady-state 
regime can directly be reached. In the present implementation, the ISPI scheme 
converges for sufficiently high temperature or sufficiently large bias voltage. The non- 
local time correlations induced by integrating out the fermions in the leads within this 
approach become long-ranged only when T and V both approach zero. Such a slow decay 
of time correlations necessitates a long memory time, i.e., one has to take into account 
more and more memory slices within the path sum in ( ITUl) . where the computational 
effort scales exponentially in the number of memory slices. This problem becomes more 
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severe for increasing U, and an exploration of the deep Kondo limit remains challenging 
at present, since the long-range correlations in time are essential for taking into account 
characteristic features of the Kondo limit. Further improvements of the scheme are 
under study to allow more efficient calculations in the large-?/ limit. However, as we 
have shown above, the present ISPI implementation is able to provide highly accurate 
and reliable results for U/T < 3, and the zero temperature limit can be reached for 
voltages V > O.IT. Moreover ISPI allows for accessing the full real-time characteristics 
of observables of interest up to in principle asymptotically long times in contrast to real- 
time QMC approaches [45 J where the error bars increase due to the sign problem for long 
times. It should also be stressed that ISPI is particularly well suited to treat the limit 
of intermediate-to-high voltages, where alternative methods often run into difficulties. 

To summarize, we presented a detailed comparison of theoretical methods designed 
to model nonequilibrium transport, taking the example of the single- impurity Anderson 
model. A very encouraging quantitative agreement between the fRG, ISPI, and tDMRG 
could be established, suggesting that a well-equipped toolbox is available to study 
nonequilibrium transport in nanostructures. We can also conclude that our methods 
fair well against alternative developments, namely TD-NRG [59J and two QMC variants 
[4"Tj SU [51] , of course keeping in mind that not all schemes can be applied in all parameter 
regimes. In conclusion, each separate approach to nonequilibrium quantum transport 
bears its own right since different aspects of the real-time dynamics of correlated 
nonequilibrium quantum systems can be learned from these different approaches. 
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Figure 1. Real-time evolution of the current I(t) at the symmetric point Eq = — U/2 
for U/T — 3 for the tDMRG. Real-time results from the ISPI approach can be found 
in [71], figure 5 therein. 
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Figure 2. Steady-state current / at the symmetric point Eq = —U/2 for U/T = 2 as 
a function of the applied bias voltage eV/T. The tDMRG results are taken from [50] 
and were computed for U = W/8 (T chosen correspondingly) and extrapolated in the 
inverse system size. The real-time QMC data are taken from [17] (U/T = 2) and [35] 
(U/T = 4 and 8). In the ISPI approach, we set the temperature to T = 0.1T. 
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Figure 3. Current-voltage characteristics (top) and differential conductance (bottom) 
for U/F = 3. The temperature is T = T/25.6 and the remaining parameters are the 
same as in figure O The circles represent the results from a QMC method using an 
analytic continuation to imaginary Matsubara voltages; [ST] (figure 2 therein). 
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Figure 4. Comparison of the steady-state current / divided by eV for U/T = 5 for 
the fRG, tDMRG and the TD-NRG 59 . The remaining parameters are the same as 
in figure [2] 
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Figure 5. Steady-state current for the mixed valence regime for different single 
electron energies E with U/T — 2. The remaining parameters are the same as in 
figured The tDMRG data are taken from [50] . 




Figure 6. Steady-state current I for finite magnetic field B, remaining parameters 
are the same as in figure [2] The solid lines are the fRG results corresponding to the 
symbols of the ISPI results and the corresponding tDMRG results are marked with 
the black stars. 




Figure 7. Temperature dependence of the steady-state current. The solid lines 
correspond to the fRG results, the symbols denote the ISPI results. The remaining 
parameters are the same as in figure [21 with j3 = X/T being the inverse temperature. 



